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The parity of a bit string of length A'' is a global quantity that can be efhciently compute using a 
global counter in 0{N) time. But is it possible to find the parity using cellular automata with a set 
of local rule tables without using any global counter? Here, we report a way to solve this problem 
using a number of r = 1 binary, uniform, parallel and deterministic cellular automata applied in 
succession for a total of 0(A''^) time. 
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PACS numbers: 89.75.-k, 05.45.-a, 05.50.+q, 89.20.Ff 



Cellular automaton (CA) is a simple local interaction 
model used to study the evolution and self-organization 
of various physical and biological systems [Q. And at 
the same time, CA can also be viewed as a restrictive 
model of parallel computation without common global 
memory. In fact, recently there is an increasing interest 
in using CA to perform certain computational tasks 
It is therefore natural to ask if CA can be used to per- 
form a task that depends on the global information of an 
input state. 

An example of this kind is called the density classifi- 
cation problem (DCP). In this problem, we are given an 
one dimensional array of bit string in periodic boundary 
conditions. We are required to apply some CA rules so 
as to evolve the state to all zero if the number of zeros 
in the input state is greater than that of ones. Similarly, 
we have to evolve the state to all one if the number of 
zeros in the input state is less than that of ones. 

After a number of fail attempts, Land and Belew 
proved that no single CA can solve the DCP without 
making some mistake And yet later on, Capcar- 
rere et al. argued that a single CA rule can solve the 
DCP provided that we modify both the required out- 
put of the automaton and the boundary conditions used 
Nonetheless, we do not completely agree with their 
approach for we have to scan through the states of the 
entire final bit string, in general, before knowing the an- 
swer. This requires either global memories or a table 
that scales with the system size in the read out; hence it 
somehow defeats the purpose of restricting ourselves to 
the use of CA in the first place. (In contrast, the read 
out in the original DCP can be determined by looking at 
any one or two bits in the final string.) Although DCP 
cannot be solved using one CA rule, Fuks showed that 
this can be done by applying two CA rules in succession. 
More precisely, he found a solution to the DCP by apply- 
ing the CA rule 184 a fixed number of times depending 
only on the string length and then followed by the CA 



rule 232 a fixed number of times again depending only on 
the string length A number of related problems de- 
pending on the global density of a string have also found 
to have CA solutions 0. 

Another challenge for CA is the parity problem (PP), 
namely, to evolve a given input bit string a using a se- 
quence of CA rules to all Po-s where is the parity of 
a. (Then the parity of the input bit string can be de- 
termined by looking at any one of the final bit string.) 
This problem appears to be much harder than the DCP 
because the output is altered simply by a flip in any one 
of the input bits. In fact, Sipper proved that no single 
r — 1 CA rule can solve the PP with fixed boundary 
conditions ||]. 

In this letter, we show that PP can be solved by apply- 
ing a number of r = 1 CA rules in succession. The term 
CA in this letter shall mean a local synchronous uniform 
deterministic binary CA rule with parallel update in pe- 
riodic boundary conditions. That is, the state of each bit 
in the next time step depends deterministically only on 
the state of its finite neighborhood. The states of all sites 
are updated in parallel and the rule table of the CA is 
covariant under translation along the bit string. Besides, 
we restrict ourselves to consider only those sequence of 
CAs that solves the PP exactly without making any mis- 
classification. 

To have a feeling of the difficulty of this problem, two 
remarks are in place. First, except for the boundary con- 
ditions, relaxing any of the above conditions makes the 
PP trivial. For example, applying the CA rule 60 first to 
the 2nd bit, then the 3rd bit and so on till the last bit 
in the bit string, the value of the last bit is the parity of 
our input string. Second, one set of local CA rules is not 
sufficient to solve the PP for input string of any length. 
The reason is simple: if such a set could solve the PP, it 
would evolve an odd parity bit string a to all one. The 
uniformity condition implies that this set would evolve 
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the concatenated even parity bit string ao to all one as 
well. But this is absurd. 

Thus, we need to invoke more than one CA rules to 
solve the PP. In fact, all we need is a few r = 1 CA rules 
reported below in Wolfram's notations 

• i?222 is the Wolfram elementary CA rule 222. 
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000 001 010 oil 100 101 110 111 



222- 



1 1 



1 



1 1 



1 



This rule replaces the two ending zeros of a string of 
zeros by ones, if the number of zeros in the string is 
more than two. Since we are using periodic bound- 
ary conditions, the (global) parity of the configura- 
tion does not change although locally parity does 
change. If we apply the rule [Af/2J times, then 
there will be no more consecutive zeros. 

^132 is similar to -R222, but it replaces the ending 
ones by zeros. It also conserves parity. 



R 



132- 



000 001 010 oil 100 101 110 111 

"0 i i~ 



• i?76 will flip a configuration of all ones to all zeros. 
000 001 010 Oil 100 101 110 111 



i?76: 



1 1 1 



• -R254 preserves the configuration of all zeros. In 
fact, careful application of i?254 together with i?76 
will put our result to the desired form. 



R 



254- 



000 001 010 oil 100 101 110 111 

"0 1 i 1 i i 1 ~ 



• i?i84 is the so-called traffic rule. It tries to move an 
one to the right if the site at right is a zero. 



R 



184- 



000 001 010 oil 100 101 110 111 
1 1 1 1 



• ^1^252 is an auxiliary rule to change a zero right of 
an one to one. 



R 



252- 



000 001 010 oil 100 101 110 111 

~o i 1 1 1 1 V 



Let a be an arbitrary input bit string of length TV, and 
Per its parity. We write Ra the resulting configuration 
after the rule R is applied to a once, for any CA rule R. 



Theorem. Let 



ai = {Ri32'--'R 



\ N \ \ N \ 



222 



Ifl 



(1) 



which means that we apply -R222 to a [N/2\ times, then 
apply i?i32 [N/2\ times and back to i?222 and so on. The 



total number of time steps is 2[7V/2J^. 
(a) If N is odd, then 



0^ if Pa = 0, 
" 1^ if Pa = 1. 



(2) 



where O''^ denotes a string of TV consecutive zeros, that 
is, the state at all the sites is zero; and similarly for 1^. 
Total number of time steps is (iV — 1)^/2. 
(b) If iV = 2q, where q is odd, then 

^ / 0~ or 1^ if Pa - 0, 

'^^ \ other than 0^ and 1^ if = 1. ^ ' 

Applying the C A rules S = P254 '^^^^^ Rje , we can trans- 
form a I to the form given in (a). 



S(7i = <^ TV 



0^ if Pa = 0, 
1^ if = 1. 



(4) 



The total number of time steps is {N'^/2 + N/2 + 1). 
(c) If = 2'" (7, where m > 2 and q is odd, let 



and 



then 



I — I I — I 

T = i?i32 ^ ^222 ^ -R184P252 



(T2 = J (Tl , 



(5) 



(6) 



^ ^ 0^ or 1^ if Pa = 0, 

^ other than 0^ and 1^ if = 1. ^ ' 



Similar to (b), S a2 is either all zeros or all ones, depend- 
ing on whether Pa is even or odd. The total number of 
steps is A^V2 +{m- 1){N + 2) + N/2 + 1. 

Remark. If we are allowed to change the lattice size, then 
we do not need (b) and (c) of the Theorem to solve the 
parity classification problem. Suppose we are given a bit 
string with even length. To find out the parity of the 
number of ones in it, we just concatenate a single zero 
bit to the bit string. The resultant bit string will be in 
odd length and hence we are back to (a) of the Theorem. 

The following two Lemmas are required to prove this 
Theorem. 

Lemma 1. For any N , ai could only be one of the follow- 
ing three forms: 0^, 1^ or (lO^'"^)*"', where l,k >1 are 
integers, 2lk = N and Pa = k (mod 2). (The notation 
^^Q2i-i^k Y<aea.ns, for example ii I — 2 and fc = 3, the bit 
string 100010001000.) 

Proof. It is obvious that both 0^ and 1^ are fixed points 
of -R222, Pi32- Thus, in the rest of this proof, we shall 
only consider configurations a with both zeros and ones. 
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As we have discussed above, for any configuration a, 
I — I 

the string R222 ^ tr has no consecutive zeros. Thus, its 
general form will look like 



R 



222 



L 2 J /T — 



a = •••101"i01"^01-' 



(8) 



where the numbers of ones between two zeros are 

greater than or equal one. (The trivial configuration 1^ 

is also possible.) 

Since i?i32 replaces the ending ones by zeros, if any 

of the rii is even, the corresponding string of ones will 

I — 

be completely annihilated after applying i?i32 ^ J . For 
example, 

i?i32^(- • • 0101111010 • • •) = i?i32(- • • 0100110010 • • •) 

= ••• 0100000010 •• • . (9) 

This example also shows that after a single pass of 
I — I I — I 

R = i?i32 2 i?222 ^ , we cannot conclude that the num- 
ber of zeros between two ones is odd. However, we can 
conclude that if there is a string of even number of zeros 
between two ones, after a single pass of R, the number of 
zeros in that string will increase by at least two. Thus, 
i?L"2"J(T must be of the form 



. . . io^"'i+^10^"^+ii . . . 



(10) 



where n[ are non-negative integers (or the trivial cases 
0^ or 1^). 

We still need to show that the n'^ are all equal. This can 
be illustrated by an example, consider the case N = 8, 



I JV I \ N \ 

Ri32 R222 (10100000) 

i?132^*J (10111011) 

(00010001) . 



(11) 



We see that the numbers of zeros between the ones tend 
to equalize, from one and five zeros to three zeros in this 
example. More precisely, R222 L"^-l (102"i+il02"2+ii . . .) 
equals 01"i+"2+i01"2+"3+io • • • up to a translation. So, 
after applying R for sufficiently long time, the resultant 
string will be in the form 1^, 0^ or Eq. (10). Since the 



first two cases are fixed points of i?, we only need to fo- 
cus on the third case. If the resultant string is in the 
third form, then it is easy to check that the total num- 
ber of ones does not increase after each application of R. 
Hence, eventually the number of ones in the string will 
stay constant under repeated application of R. And this 
happens if and only if n[ are all even or all odd. In this 
case, 

i?(102"'i+il02"^+H---) 

up to a translation. Hence, the repeated application of 
R equalizes the number of zeros between the ones; and a 



configuration in the form (10^"+^^)' is a fixed point of R. 
Finally, it is straight forward to check that at most [-^J 
applications of R is enough to bring a bit string to the 
fixed points of the form 0^, 1^ or (lO^'-i)'^. (One of the 
worst cases is the configuration 1(10)^^^^^^^ for N odd.) 

Since i?i32 and i?222 conserve parity, k (mod 2) is 
equal to the parity of (10^'~^)''', which is just P^- This 
completes the proof of Lemma 1. □ 

Lemma 2. With the notations of the Theorem, we have 



rp _ 



T(102'-i) 



(lO'-i) 



2k 







I A' 

N 



if I is even, 

if ? = 1, 

if Z > 3 is odd, 



(13) 
(14) 

(15) 



where fc > 1. 

Proof. Eq. (^ and Eq. (|l|) are trivial. For Eq. (|l|), we 

consider two cases, lil > 2, 

T(102'-l)'== fi?132L^Ji?222^*Ji?184i?252) (10^'-^)'= 



= (i?132L-Ji?222L-Ji?184) (llO^'-^)'^ 
'i?132L^Ji?222L*j) mO^'-')' 



(lO'-ilO'-i)*-' if Z is even, 
(02')'= if I is odd. 



(16) 



If ; = 1, 



T(IO)''^ = (i?132L^Ji?222L*Ji?184i?252) (10)* 
= fi?132L^Ji?222L^Ji?184l l"^ 



= 1 



N 



This concludes the proof of Lemma 2. 



(17) 

□ 



Proof of the Theorem, (a) If N is odd, among the 
three forms provided by Lemma 1, (lO^'"^)*^ cannot be 
reached. Since parity is conserved, the final configuration 
can only be 0''^ or l'^ according to the initial parity. 

(b) If A^ = 2q, parity of O'^ and l'^ are even, but the 
parity of (lO^'"^)*^ must be odd because k is odd. Thus, 
we have proved Eq. (||). Now, by flipping the all ones to 
all zeros and changing all other configurations to all ones 
by the rules 5", we put the final configuration to a form 
similar to (a). 

(c) If P„ = 1, then = 1 A: (mod 2) and 
/ = 2™^^(7' where q' is odd. By Lemma 2, we have 

/j^m — 1 ^-j^Q2/ — rpm—l^-^Q2^q' — l'jk 



(10V-i)2 



(18) 
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If Per = 0, then I = 2™ ^q' where m' < m and q' is 
odd. The above equation does not reach the last hne and 
rpm-i(^lQ2i-i-^k jg qN qj. ;^jv rpijig coHipletes the proof of 

the Theorem. □ 
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To summarize, we show that it is impossible for any 
single set of CA rules to correctly compute the parity of 
a bit string. But surprisingly, we find a CA solution to 
the PP using a sequence of r = 1 CAs. For an input 
bit string of length N, the worst case run time scales as 
0(-/V^). However, we have no idea if the present method 
has the shortest worst case run time or not. And since 
the PP is equivalent to the computation of the sum of a 
given bit string modulo two, our result implies that CA 
can be used to count the number of ones in a bit string 
modulo two. In other words, if we denote the number of 
ones in a bit string a by ^^cr, then one can use a sequence 
of CAs to compute the least significant bit of #(7 written 
in binary notation. Besides, the CA solution of the DCP 
for a bit string 2" — 1 long can then be regarded as a 
way to compute the most significant bit of #cr. There- 
fore, it is instructive to investigate the possibility of using 
CAs to compute; any given bit of =ff(7 and hence to count 
the number of ones in a bit string provided that we have 
0(logA'') copies. Along a similar line of thought, it is 
also worthwhile to look for CA solutions to the problem 
of addition over a finite ring or field. 
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